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Определение интервалов квазистационарности 
экономических систем 


В работе рассмотрен вопрос определения оптимального интервала адаптации алгоритма динамического 
управления капиталом для нестационарного случая методами расчета показателя Херста и построения 
автокорреляционной функции для анализа временных рядов. Проведен анализ влияния выбора интервала 
адаптации на эффективность алгоритма. Из анализа полученных результатов следует, что метод расчета 
показателя Херста позволяет более эффективно, чем метод построения автокорреляционной функции, 
определить интервал стационарности модели функционирования экономической системы. 


Введение 


В своей монографии [1] Р. Винс разработал оригинальную теорию «оптималь- 
ного }›. В работах [2-4] нами показано, что применение эмпирического метода нахожде- 
ния «оптимального }» Р. Винса в алгоритмах динамического управления капиталом 
(ДУК) нуждается в серьезной корректировке. Применение новых методов позволяет 
существенно повысить эффективность алгоритмов ДУК. Это наглядно проявляется 
для нестационарных моделей функционирования экономических систем. В то же 
время остается открытым вопрос о влиянии интервала адаптации на эффективность 
предлагаемых алгоритмов ДУК и выбора оптимальной величины этого интервала. 


Анализируемые модели функционирования 
экономической системы 


Целью исследований в данной работе является анализ эффективности алго- 
ритма ДУК при выборе различных интервалов адаптации. Эффективность алгоритма 
оценивалась величиной среднего геометрического (С’) доходности, достигнутой для 
алгоритма за отчетный период. Задачей оптимизации являлось достижение макси- 
мального значения величины С’для анализируемой модели функционирования эконо- 
мической системы. Одной из задач исследования была разработка метода определения 
интервала квазистационарности методом анализа показателя Херста временного ряда. 

В качестве экономико-математической модели функционирования экономичес- 
кой системы, как и в работах [2-4], используется модель, входными характеристи- 
ками которой является одномерный временной ряд случайных величин выигрышей и 
проигрышей ( {х,} ). 

Пусть х, — значение 1-го элемента функционирования экономической системы. 


Предположим, что на интервале стационарности временной ряд этих значений мо- 
жет быть описан законом распределения с функцией плотности вероятностей И’(х) 


(здесь х — значение случайной величины). 
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В большинстве экономических систем значения {х‚} подчиняются нормальному 
закону распределения, что вытекает из центральной предельной теоремы. Таким обра- 
зом, на интервале стационарности модель может быть задана двумя числовыми харак- 
теристиками: математическим ожиданием (М) и дисперсией (52°). В то же время 
для реальной экономической системы на значения {х,} накладываются некоторые 
ограничения, которые в нашем случае формулируются следующим образом: 

— максимальный проигрыш не может быть больше 100% (таким образом, х; > -1\1); 
— максимальный выигрыш не может превышать некоторого значения (в работе при- 
нято х, <1 У). 

Для генерации данных использовались стандартные средства табличного процес- 
сора М$ Ехсе|, алгоритм подробно описан в работах [2-4]. Для имитации нестацио- 
нарности реализации {х‚} проводилась модуляция ряда по синусоидальному закону с 
частотой © («медленные» изменения), 100 («быстрые» изменения) и без модуля- 
ции (стационарный ряд). При генерации модулированных данных использовались 
два различных интервала стационарности: т„ = 20 («короткий» интервал) и т.„ = 60 
(«длинный» интервал). Таким образом, для характеристики нестационарной модели 
функционирования экономической системы наряду с параметрами М и 5Р вводятся 
параметры т„ и О. Были исследованы 5 различных моделей функционирования 
экономической системы: 1) стационарная модель (СМ); 2) модель с «медленными» 
изменениями и «длинным» интервалом стационарности (МДМ); 3) модель с «мед- 
ленными» изменениями и «коротким» интервалом стационарности (МКМ); 4) мо- 
дель с «быстрыми» изменениями и «длинным» интервалом стационарности (БДМ); 
5) модель с «быстрыми» изменениями и «коротким» интервалом стационарности (БКМ). 


Анализируемые алгоритмы ДУК 


Вопросы выбора алгоритма ДУК подробно обсуждались нами в [2-4]. Данная 
работа посвящена анализу влияния выбора интервала адаптации на эффективность ал- 
горитма, в связи с чем, вопрос о выборе анализируемого алгоритма не является акту- 
альным. Для анализа был выбран эмпирический метод «оптимального /» Р. Винса [1], 


который может быть описан формулой: 


Кия Ру" | (1) 


= Х тт 


где п=т.„_ интервал адаптации для оценки последующего поведения ряда, Хх — ве- 
личина выигрыша или проигрыша системы (исторические данные), взятая с про- 


тивоположным знаком; хи — самый большой проигрыш системы (всегда со знаком 
«—), вычисленный на интервале адаптации. 

Метод определения Г” для каждого шага подробно описан в [3] и сводится к 
так называемому «наивному прогнозу» [5]. В этом случае считается, что если на п-м 
шаге адаптации определена величина Г °, то на (п+1)-м шаге она изменяется незна- 
чительно и ее можно считать прогнозной. 

Вопросу выбора т.„ внимания не уделялось, использовалось известное значе- 


ние интервала стационарности анализируемой модели функционирования экономи- 
ческой системы. Следует отметить, что определение величины этого интервала не 
является тривиальной задачей. 
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В то же время можно предполагать, что выбор малого интервала адаптации 
<^ * 
алгоритма и <<т,„ приводит к невысокой точности прогнозного значения /, а 


выбор слишком большого интервала адаптации п>> т. приводит к искажению 


* 
прогнозного значения Г. 


При проведении оценки алгоритма ДУК в качестве критерия эффективности 
была использована величина С’, достигнутая для алгоритма за отчетный период ( [): 


б=(К/Ко)"", (2) 


где Ки А — соответственно величина наращенного капитала в конце отчетного пе- 
риода и величина начального капитала инвестора, Г — количество сделок. 


Результаты анализа моделей функционирования 
экономических систем методом построения 
автокорреляционной функции 


Стандартным методом анализа динамических рядов является метод построения 
автокоррелограмм [6]. Автокоррелограмма показывает численно и графически авто- 
корреляционную функцию (АКФ) — коэффициенты автокорреляции для последова- 
тельности лагов из определенного диапазона. Интерес в этом случае представляют 
наиболее сильные (высоко значимые) автокорреляции. 

Анализ АКФ для МКМ системы свидетельствует о высоком значении показа- 
теля автокорреляции для лагов 12 <т<6| и 93<т1<180 (отличен от 0 на уровне 
значимости р< 0,05). Для МДМ системы высокое значение показателя автокорре- 


ляции для лагов |< т < 140. Для СМ системы высокое значение показателя автокор- 
реляции (отличен от 0 при р < 0,05) наблюдается для лагов 126 <т<178. 

Таким образом, можно полагать, что для МКМ системы интервал адаптации не 
должен превышать т = 60, а для МДМ систем интервал адаптации не должен пре- 
вышать 7 = 140. 

Анализ АКФ функции для БКМ системы не позволяет сделать какого-либо за- 
ключения о характере динамического ряда. Анализ АКФ для БДМ системы свиде- 
тельствует о высоком значении показателя автокорреляции для лагов 1<1<80 и 
140 <г< 200 (отличен от 0 при р < 0,05). 


Алгоритм анализа нестационарных временных рядов 
с помощью показателя Херста 


Для характеристики поведения нестационарных временных рядов в настоящее вре- 
мя применяется метод расчета фрактальной размерности ряда (показатель Херста Н) 
[7 [8]. Значение показателя Херста может изменяться от 0 до 1. При Н =1 ряд 
называется персистентным (последующие изменения значений временного ряда 
наследуют его предыдущее поведение), при Н = 0,5 ряд будет случайным (после- 
дующие значения временного ряда не связаны с его предыдущими значениями), при 
Н = 0 ряд будет антиперсистентным (последующие изменения значений временного 
ряда противоположны его предыдущему поведению). При этом характер поведения 
ряда может меняться на различных временных промежутках. 
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В работе [9] представлен один из методов расчета пока- 
зателя Херста временного ряда, известный как Оефепаеа ЕТас- 
бпаноп Апа1уз15 (ОЕА). Данный метод позволяет рассчитать 
показатель Херста для различной глубины анализа т. 

При проведении расчета методом ОЕРА весь временной 
ряд анализируемых показателей Х’ разбивается на отрезки рав- 
ной длины т. Для каждого К -го отрезка разбиения определя- 
ется линейный локальный тренд 2(1) =а-Е+Ь. Для этого от- 
резка ОЕА функция Р определяется соотношением: 

(+1 
Ера)=- > @-20). (3) 
1=КТ+1 

Среднее значение этой функции < Е°(т)> по всем [/т 


интервалам разбиения есть функция от т: 


а 


1 
Ввод нестационарной 
реализации {х;} 


Начало процедуры 
анализа, =4 


3 


Выбор интервала 
разбиения = 


4 
Вычисление 


<Е* (2) по 3) и 


6 
Расчет показателя 


г = 
ВЫ АЯ <Е*(г)>=-У ЕР (Ф). (4) 
Г К=0 
При этом < Ё*(т) > связано с г соотношением [8]: 
Ви, (5) 
Рисунок 1 — где Н - показатель Херста ряда. Отсюда, строя зависимость 


Структурная схема 
реализации метода 
ОКА расчета щеи прямои и оценивается значение показателя Херста по алго- 


показателя Херста ритму, приведенному на рис. 1. 


№2( <Е”(г)> ) от (г), определяется наклон аппроксимирую- 


Результаты анализа моделей функционирования 
экономических систем методом РЕА 


На рис. 2 в качестве примера представлены результаты анализа фрактальной 
размерности МДМ модели функционирования экономической системы. 

Из проведенного анализа следует, что для МДМ системы для интервалов 
т<48 характер поведения ряда близок к случайному (Н = 0,5), для интервалов 
48 <г< 209 ряд является антиперсистентным (Н =0), для интервалов т > 209 ряд 
является персистентным (Н =1). Таким образом, для МДМ системы интервал 
стационарности 48 <т„ < 209 (значение при генерации т„„ = 60). 

Для МКМ системы для интервалов т <22 характер поведения ряда близок к 
случайному (Н = 0,5), для интервалов 22 <г< 63 ряд является антиперсистентным 
(Н =0), для интервалов т > 63 ряд является персистентным (Н =1). Антиперсис- 
тентность ряда свидетельствует о возврате его значений обратно при их отклонении 
от среднего, что является признаком его стационарности, следовательно, интервал 
стационарности этой системы 22 <т„„ < 63 (значение при генерации т.„ = 20). 

Для СМ системы интервал антиперсистентности выделен быть не может, для 
всех значений т характер поведения ряда близок к случайному (Н = 0,5). 

В случае моделей функционирования экономических систем с «быстрыми» изме- 
нениями интервал антиперсистентности короче и выделяется не так явно, как для слу- 
чая «медленных» изменений. Для БКМ системы интервал стационарности 24 <т„<29 


(значение при генерации т, =20), для БДМ системы интервал стационарности 
43 <г.„ <110 (значение при генерации т.„ = 60). 
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ш(<Е? (т) 1?) 


'л 


у = 1.35х- 6.72 


у =0.60х - 1.79 ы(5) 


Рисунок 2 — Результаты анализа фрактальной размерности временных рядов 
МДМ модели функционирования экономической системы методом ОКА. 


Сравнивая полученные результаты с результатами для метода анализа АКФ, 
можно отметить большую эффективность метода расчета фрактальной размерности 
ряда, который позволяет оценить верхнюю и нижнюю границу интервала квазиста- 
ционарности ряда. 


Результаты анализа алгоритма ДУК для различных 
интервалов адаптации 


Для проверки предположения о влиянии интервала адаптации алгоритма ДУК 
на его эффективность были проведены расчеты значений С’„ для различных 
значений (ип) этого интервала. На рис. 3 в качестве примера приведена зависимость 
С „(п) для МДМ модели функционирования экономической системы. 


| С опт. > 
1.19 - 


1.18 - с 


1.17 


0 40 80 120 160 200 240 280 
п 


Рисунок 3 — Зависимость С’, от значений интервала адаптации и 


опт 


для МДМ модели функционирования экономической системы 


Из анализа полученных результатов следует, что первоначальное предположе- 
ние о наличии оптимального интервала адаптации п, которое не должно сильно 
отличаться от интервала стационарности т,„, верно. Так, для МКМ системы макси- 


ст? 


мального значения С’„ достигает при п = 80, для МДМ системы -— при п = 140. 


опт 
Аналогичные результаты получены и для моделей функционирования экономи- 
ческих систем с быстрыми изменениями. 
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Следует отметить, что анализ, проведенный для алгоритмов более эффектив- 
ных, чем эмпирический метод определения «оптимального /» Р. Винса, дает более 


выраженную зависимость С’ от и. 


опт 


Выводы 


Из проведенного анализа можно сделать следующие выводы: 
1. Метод построения АКФ позволяет дать верхнюю границу оценки интервала 
квазистационарности т.„„, модели функционирования экономической системы. 


ст ? 


2. Метод ОРА анализа временных рядов с помощью показателя Херста позво- 
ляет эффективно определить интервал стационарности т„ модели функционирова- 


ния экономической системы и является в этом отношении более эффективным, чем 
метод построения АКФ. 

3. Эффективность метода ОКА анализа временных рядов с помощью показате- 
ля Херста снижается при повышении степени нестационарности временных рядов, 
оставаясь при этом более эффективным, чем метод построения АКФ. 

4. Эффективность алгоритма ДУК зависит от интервала адаптации и и близка 
к оптимуму при значениях и, близких к т». 
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Т.В. Гур’янова 

Визначення нтервалв квазстацюнарност! економ!чних систем 

Робота присвячена питанню визначення оптимального 1нтервалу адаптацй алгоритму динам!чного 
керування кашталом для нестацонарного випадку за допомогою метод1в розрахунку показника Херста 1 
побудови автокорелящшино! функцй задля аналйзу часових рядв. Проведено анал!з впливу вибору 
1нтервалу адаптаци на ефективнисть алгоритму. Порвняння результатив проведеного анал!зу дозволяе 
стверджувати, що метод розрахунку показника Херста дозволяе блыш ефективно, иж метод побудови 
автокорелящйно! функций, визначити 1нтервал стацюнарност! модел! функцюнування економ1чно] системи. 


Статья поступила в редакцию 28.12.2009. 
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